Magnetic moment of an electron near a surface with dispersion 
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Boundary-dependent radiative corrections that modify the magnetic moment of an electron near 
a dielectric or conducting surface are investigated. Normal-mode quantization of the electromag- 
netic field and perturbation theory applied to the Dirac equation for a charged particle in a weak 
magnetic field yields a general formula for the magnetic moment correction in terms of any choice 
of electromagnetic mode functions. For two particular models, a non-dispersive dielectric and an 
undamped plasma, it is shown that, by using contour integration techniques over a complex wave 
vector, this can be simplified to a formula featuring just integrals over TE and TM reflection coeffi- 
cients of the surface. Analysing the magnetic moment correction for several models of surfaces, we 
obtain markedly different results from the previously considered simplistic 'perfect reflector' model, 
which is due to the inclusion of physically important features of the surface like evanescent field 
modes and dispersion in the material. Remarkably, for a general dispersive dielectric surface, the 
magnetic moment correction of an electron nearby has a peak whose position and height can be 
tuned by choice of material parameters. 



The spin magnetic moment of an electron differs from 
the value predicted by the Dirac equation because the 
electron is coupled to the quantized electromagnetic field. 
The electromagnetic field is in turn modified by the pres- 
ence of a material boundary; thus the magnetic moment 
for an electron near a surface differs from the value found 
in free space (J-Q. This shift is, in principle, measur- 
able spectroscopically, however previous investigations 
have shown that this effect is so minute as to be im- 
measurable with currently attainable experimental pre- 
cision. All these previous works [U-Q model the surface 
as a perfect reflector, which forces the parallel electric 
and the perpendicular magnetic fields on the surface to 
vanish and can be represented by a medium with infinite 
refractive index. While simple to work with, this model 
has some obvious physical deficiencies in that it ignores 
electromagnetic field modes that are evanescent outside 
the medium and in that the medium does not become 
transparent at infinite frequencies as any real material 
would. For some quantum electrodynamic surface effects 
like the Casmir-Poldcr force on an atom in front of a sur- 
face Q, the idealized model of a perfect reflector works 
just fine and its results are reproduced as limiting cases 
e.g. of a non-dispersive dielectric medium but with in- 
finite refractive index 0] or a dispersive plasma surface 
but with infinite plasma frequency uj p (lo| . This Letter is 
going to show that the situation is radically different for 
the spin magnetic moment of an electron near a surface 
where different models for the electromagnetic response 
of the surface give quite different results, not necessarily 
obtainable as limiting cases of each other. That this can 
happen was already found for the inertial mass shift of 
a free electron near a surface 111, the reason for which 



was found to be connected to the fact that a free parti- 
cle admits excitations of arbitrarily low frequency, thus 
interacting with photons of arbitrarily long wavelengths 
which arc dealt with very differently in dielectric and 



perfect-reflector models of a reflecting surface. 

In order to determine the quantum-field theoretical 
corrections to the electron's magnetic moment, one usu- 
ally calculates the vertex diagram to the respective order 
of interest. To 1-loop order e 2 = a and in free space, 
this is a straightforward calculation and graduate text- 
book material (c.f. e.g. [l2|), but when the electromag- 



netic field is reflected by a surface even 1-loop calcula- 
tions of quantum electrodynamics get rather complicated 
fl3j ] , largely because of the loss of translation invariance. 
Here we are interested only in the correction to the mag- 
netic moment due to the presence of the material surface 
nearby, whence we take a different, more appropriate ap- 
proach that dispenses with quantizing the electron field 
and just works with the quantized photon field and (first- 
quantized) electrons. This works because the Feynman 
graphs that contribute to 1-loop order are the irreducible 
vertex correction and reducible graphs of the bare ver- 
tex with 1-loop self-energy and vacuum-polarization in- 
sertions (cf. e.g. Fig. 12.15 of (l3|). and of those the 
vacuum polarization is the only contribution that would 
necessitate the second-quantization of the electron. How- 
ever, the vacuum polarization does not, to 1-loop or- 
der, contain any internal photon lines and hence to this 
order does not receive any boundary-dependent correc- 
tions. Thus, we can work out the boundary-dependent 
part of the shift by quantizing just the electromagnetic 
field and use perturbation theory in the Dirac equation 
to determine the energy shift — fi ■ Bo, due to the pres- 
ence of the material surface, of the electron's spin in a 
weak external magnetic field Bo- Alternatively, one can 
of course take a non-relativistic expansion of the Dirac 
equation first, by means of a Foldy-Wouthuysen transfor- 
mation or or otherwise, and then use perturbation theory 
to work out the same energy shift, as has been done for 
working out the magnetic moment correction of the elec- 
tron near a perfect reflector This approach requires 
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care since several successive orders in the non-relativistic 
expansion turn out to contribute to the shift Q. 

We start by considering a quantized electromagnetic 
field interacting with an electron that sits near a half- 
space of constant refractive index n. Initially we shall 
assume that the weak classical magnetic field Bo = Bqz 
interacting with the magnetic moment is normal to the 
surface spanning the xy plane, so that in the resulting en- 
ergy shift the coefficients of terms linear in a z Bo yield the 
spin magnetic moment. We write the quantized electro- 
magnetic field in terms of mode functions fkA and stan- 
dard photon creation and annihilation operators as 

A= ^ (f k Aa k A + fkA«L): (1) 

all modes 

where A denotes a particular polarization, and any nor- 
malization constants have been absorbed into the mode 
functions. We then apply perturbation theory to the 
Dirac equation for an electron coupled to the electromag- 
netic field via minimal coupling by using the electron's 
states in a constant magnetic field [HJ as unperturbed 
states and the coupling term — e^oj ■ A with the quan- 
tized photon field as perturbation. Second-order pertur- 
bation theory gives relative magnetic moment corrections 
A/z//z of order e 2 . The non-relativistic expansion we 
make in the process eventually turns out to lead to a 
series in powers of the ratio of the electron's Compton 
wavelength and its distance z from the surface, 1/mz 
[l6| . The leading order boundary-dependent correction 
to the magnetic moment is of order A/1//1 pa e 2 /(mz) , 



which in fact relates to another main reason, beside tech- 
nical complications, for using an approach based on the 
quantized photon field interacting with a first-quantized 
Dirac electron over a fully second-quantized quantum 
electrodynamic approach: as shown by Kreuzer Q for a 
perfectly reflecting surface, the quantum electrodynamic 
approach gives leading-order contributions to the g fac- 
tor as well as to the mass shift that are of order e 2 /(mz) 
but which cancel in the boundary-dependent corrections 
to the magnetic moment. The same happens if dielectric 
models of the surface are used [l7j], and thus the deriva- 
tion of the boundary-dependent magnetic-moment cor- 
rections requires the notoriously difficult calculation of 
next-to-leading-order terms in 1/mz. By contrast, per- 
turbation theory applied to the Dirac equation, as we use 
here, or to the Schrodinger equation, as used in Ref. 0], 
circumvents any such problems because one evaluates the 
magnetic moment shift directly and the leading terms are 
of order e 2 /(mz) 2 from the outset. 

The application of perturbation theory is straightfor- 
ward, if somewhat lengthy, with the one exception that 
the dipole approximation commonly used in this ap- 
proach well-known from quantum optics applications is 
not universally applicable. In some terms one needs to 
go beyond the dipole approximation by not just taking 
the value of the quantum field A at the position of the 
electron but by Taylor-expanding it around this position 
and treating the resulting powers of the displacement op- 
erator r — ro as operator acting on the Landau states of 

the electron in the external magnetic field, which gener- 

— 1/2 

ates additional factors of B ' . The end result turns 
out to be independent of the Landau level. We find 
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which, wc emphasize, is a completely general result for 
any quantized photon field of the form (TTJ). For a non- 
dispersive dielectric half-space the mode functions can be 
found in Ref. [l8|. Substituting these into the above ex- 
pression, and renormalizing the shift by subtracting the 
value it would take in free space, we find that if we write 
the photon momentum as k = k|| + k z z we can deform 



the contour in the complex k z plane to give the results in 
terms of double integrals over the reflection coefficients 
at the surface. By cycling the Cartesian coordinates in 
the mode functions we can also consider surfaces parallel 
to the external field without much additional effort and 
thus give the results for the shifts along both directions: 
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Here £ = zo; originates from the photon frequency and 
= k z /ui is the cosine of the angle of incidence. For 
later convenience, we have written the dielectric function 
e(£) as a quantity that could depend on frequency, but for 
the non-dispersive case we simply have e(£) = n 2 = e(0). 
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The terms outside the integrals come from the subtrac- 
tion and separate consideration of the point at k z = iku , 
where TM modes contribute an electrostatic term, which 
is mathematically expedient. For the non-dispersive di- 
electric these integrals can be evaluated analytically and 
yield: 
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For large refractive indices the reflection coefficients §5§ 
become those of a perfect reflector, whence one might 
expect that in the limit n — > oo the above results for A/i 
would reproduce those previously derived for an electron 
near a perfectly reflecting surface 0- However, in that 
limit we find 
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0(l/n)) (8) 



0(l/n) 



(9) 

The linear rise with n of the leading terms appears to 
be unphysical, as this would suggest that the magnetic 
moment could be increased arbitrarily by increasing the 
refractive index n of the surface. However, as we shall ex- 
plain below, this misconception arises from the unphys- 
ical assumption of a dispcrsionless medium. A curious 
observation to note is that the next-to-leading terms in- 
dependent of n in Eqs. (0) and (j9|), if taken on their own, 
do in fact reproduce the results of the perfect-reflector 
case (cf. Eq. (7.12) of )• A consistency check of tak- 
ing the limit n — > oo in the reflection coefficients (|5j| first 
and evaluating the integrals ([3j> and (j4|) afterwards repro- 
duces the results of Ref. Q and reveals the mathematical 
origin of the discrepancy: the limits n — > oo and r\ — > oo, 
which corresponds to the static limit k z — > iku , in the TE 



r 



reflection coefficient do not commute. 

To resolve the problem, we need to consider the inter- 
action of the an electron with a dispersive surface. Wc 
start with the simplest model of a dispersive medium, a 
plasma surface [19( with the dielectric function 



c( W ) 



(10) 



where u p is the plasma frequency. The normal modes AcA 
for this system can be derived following Refs. [l(], fl9j ) . 
though for the present purpose it is expedient to sepa- 
rate left- and right-incident modes and write the modes 
in the same form as for the non-dispersive medium 
On substitution of these into Eq. © we find that manip- 
ulating the contour in the complex k z plane (noting in 
particular that the contribution from a pole in the TM re- 
flection coefficient cancels with the contribution from the 
surface plasmons, as seen in [2(|) gives precisely Eqs. ([3]) 
and just with a different dielectric function entering 
into the reflection coefficients. However, for this partic- 
ular dielectric function there is a technical problem with 
the evaluation of the integrals for the TE modes. For 
this case we in fact deform the contour back down to the 
real k z axis and evaluate the integral directly. This is 
relatively straightforward because the TE reflection co- 
efficient for a plasma surface takes a particularly sim- 
ple form. The magnetic moment shifts wc obtain then 
agree in the limit u> p — > oo with the perfect-mirror re- 
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suits given by Q, but not with the limit n — > oo of the 
non-dispersive dielectric model, as explained above. This 
is not altogether surprising since neither the plasma nor 
the perfect-reflector models permit modes that are trav- 
elling inside the medium and evanescent on the vacuum 
side, which, however, feature in the dielectric model. We 
also find the at first sight worrying result that as we send 
uip — > the shift in the plasma model diverges, while one 
would expect it to be zero. Mathematically this arises 
because the u> p — > and k z — > ifcii limits of the TM re- 
flection coefficient do not commute, indicating that the 
static polarizability of the medium plays a major role in 
the shift. 

To resolve these issues we turn to a dispersive dielec- 
tric as model for the surface. Along the lines of 21 we 



introduce a restoring force into the equation of motion 
governing the dielectric response of the material, which 
leads to the dielectric function 
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(11) 



This model now includes total internal reflection inside 
the medium and evanescent modes in vacuum, so we 
expect results with features in common with the non- 
dispersive case. A derivation of the magnetic moment 
shift from first principles is now much more complicated, 
as the equation for the field modes f^x cannot be written 
as an Hermitean eigenvalue problem. Instead one would 
need to include both dispersion and absorption into the 
model by constructing a Huttner-Barnett-type field the- 
ory for the dielectric. A short-cut would be provided by 
the Lifshitz theory, where the analogy of atomic energy 
level shifts near surfaces [22[ suggests that one gets the 
same expressions © and (jH), as before, for the magnetic 
moment shifts. Substituting the dielectric function (jlip 
into Eqs. ([3H5|), we find that for a dispersive dielectric the 
magnetic moment shift exhibits a peak at certain values 
of the parameters. In Fig. [1] we plot the magnetic mo- 
ment shift near a dispersive dielectric as a function of 
the static dielectric susceptibility x(0) = e(0) — 1, which 
equals Wp/wj for a dispersive dielectric and n 2 — 1 for the 
non-dispersive model. If we were to continue the plot 
to large values of x(0), the graphs for the two models 
converge into one linearly rising graph. The peak is well- 
defined at values of \u>tz\ < 0.07 and \u>tz\ < 0.25 for Bo 
perpendicular and parallel to the surface, respectively. 
The ratio of the height of the peak to the non-dispersive 
result at the same x(0) is given by 



A/zj_disp 30.3 eVnm 



AM||dis F 



A/^_Lnondisp 



\u>TZ\ 



Ami 



nondisp 



81.6 eVnm 

\u>tz\ 



(12) 

As \u)tz\ decreases, we find that the peak moves closer 
to \/x(0) = lo p /ut ~ 2. Typical values of lot for metals 
are on the order of a few cV (sec, for example, (23|), 
meaning that an extremely small z would be required to 
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FIG. 1. A/ix for a dispersive surface with \ojtz\ = 0.02 and 
for a non-dispersive surface 



get a significant enhancement. However, materials can 
be tuned to have lot significantly smaller than this - for 
example an InSb semiconductor grating can have lot (and 
LOp) in the range of a few meV [24j. For z of a few tens 
of nanometres this corresponds to an enhancement factor 
on the order of 10 3 . 

The large x(0) behaviour of the result sheds light on 
the apparent problem of the non-dispersive result in the 
limit n — > oo. For large x(0), we observe that the shift for 
the dispersive dielectric model becomes linear in \/x(0) 
and agrees with the non-dispersive results; so for large 
x(0) the two models are equivalent. However, an arbi- 
trarily large x(0) would model a dielectric with an arbi- 
trarily large static polarizability which is physically not 
possible. The comparison with the plasma and perfect 
reflector models shows that the physics is different for 
those: in the plasma model the static limit is infinite 
right from the start, whence results do agree with the 
perfect reflector in the limit lo p — > oo. 

Mathematically, the discrepancies between the various 
results come from non-commutation of certain limits of 
the reflection coefficient, namely a non-commutation be- 
tween the static limit (fc z — > ik\\) and whichever limit 
we have to take in order to compare models. We have 
already discussed the n — > oo limit of the non-dispersive 
model. Similarly, the limit lot (x(0) — > oo) in the 
dispersive dielectric model does not reproduce the results 
of the plasma model because the limits lot — > and 
k z — > iku of the dispersive TE reflection coefficient do 
not commute. 

The shift for the non-dispersive dielectric model has a 
distance dependence of 1/z 2 for all n. However, in the 
plasma model we find a distance dependence of 1/z 3 for 
LOpZ <C 1, namely: 



Afi±(\w p z\ < 1) 



A/i||(|w p z| < 1) 



64itV2lj p 



5e J 



128ttV2w 



(13) 



(14) 
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This is because at small distances the dominant effect is 
the electrostatic interaction of the magnetic moment with 
the surface plasmon and the 1/z 3 terms come entirely 
from the surface plasmon part of the mode functions, 
similarly to what is discussed in Ref. [20| , For the same 
reasons, the dispersive dielectric model also yields a 1/z 3 
dependence of the shift in the limit lu p z <C 1. 

In summary, we have shown that the magnetic mo- 
ment shift near a dispersive dielectric surface is notably 
different from the predictions made for a perfectly re- 
flecting surface. For the case of the magnetic field Bo 
normal to surface even the sign is different than pre- 
dicted by a perfect-reflector model. Additionally, we 
have shown that the inclusion of dispersion can signif- 
icantly modify the magnitude of the effect and its dis- 
tance dependence, and that these modifications can be 
tuned by choice of material parameters. A physically re- 
alistic choice of model for the low-frequency behaviour of 
the dielectric function is essential for obtaining physically 
correct results. Our formulae dSHHJ) can be applied to 
get the magnetic moment shifts near a surface provided 
the frequency-dependence of the reflection coefficients is 
known. 

It is a pleasure to thank Robert Zietal for discussions. 
Financial support from the UK Engineering & Physical 
Sciences Research Council is gratefully acknowledged. 



[1] E. Fischbach and N. Nakagawa, |Phys. Rev. D] 30, 2356 
(1984) 



[2] K. Svozil, |Piiys. Rev. Lett. 1 54, 742 (1985) 

[3] M. BordagT TPhys. Lett. B| 171, 113 (1986) 

[4] D. G. Boulware, L. S. Brown, and T. Lee, |Phys. Rev. D| 

32, 729 (1985) 

[5] M. Kreuzer and K. Svozil, |Phys. Rev. D| 34, 1429 (1986) 
[6] M. Kreuzer, |j~Phys. A: Math. Gen.| 21, 3285 (1988) 
[7] G. Barton and N. Fawcett, Physics Reports 170, 1 (1988) 
[8] H. B. G. Casimir and D. Polder, [Phys. Rev.| 73, 360 

(1948) 

[9] S. T. Wu and C. Eberlein, Proc. R. Soc. Lond. A 455, 
2487 (1999) 

[10] M. Babiker and G. Barton, J. Phys. A: Math. Gen. 9, 
129 (1976) 

[11] C. Eberlein and D. Robaschik, |Phys. Rev. Lett.| 92, 
233602 (2004) 

[12] M. E. Peskin and D. V. Schroeder, An Introduction to 
Quantum Field Theory (Westview Press, Boulder, 1995) 

[13] C. Eberlein and D. Robaschik, |Phys. Rev. D | 73, 025009 
(2006) 

[14] M. Le Bellac, Quantum and Statistical Field Theory 

(Clarendon Press, Oxford, 1991) 
[15] M. H. Johnson and B. A. Lippmann, |Phys. RevT| 76, 828 

(1949) 

[16] We work in natural units c = 1 = H, eo = 1. 
[17] R. Zietal, C. Eberlein, and D. Robaschik, (unpublished) 
[18] C. K. Carniglia and L. Mandel, |Phys. Rev.T)| 3, 280 
(1971) 

[19] J. M. Elson and R. H. Ritchie, |Phys. Rev. 13] 4, 4129 

(1971) 

[20] M. Babiker JPhys. Rev. B|13, 305 6 (1976) 
[21] G. Barton. iProc. R. Soc. Lond. AI453, 2 461 (1997) 
[22] J. M. Wylie and J. E. Sipe, |Phys. Rev. A|30 1185 (1984) 
[23] J. P. Walter and M. L. Cohen, |Phys. Rev. B| 5, 3101 (Apr 

1972) 

[24] J. Gomez Rivas, M. Kutt ge, H. Kurz, P. H aring Boli- 
var, and J. A. Sanchez-Gil, l |Appl. Phys. Lett.| 88, 082106 
(2006) 



